Regular networks of Luttinger liquids 
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We consider arrays of Luttinger liquids, where each node is described by a unitary scattering 
matrix. In the limit of small electron-electron interaction, we study the evolution of these scattering 
matrices as the high-energy single particle states are gradually integrated out. Interestingly, we 
obtain the same renormalization group equations as those derived by Lai, Rao, and Sen, for a 
system composed of a single node coupled to several semi-infinite ID wires. The main difference 
between the single node geometry and a regular lattice is that in the latter case, the single particle 
spectrum is organized into periodic energy bands, so that the renormalization procedure has to 
stop when the last totally occupied band has been eliminated. We therefore predict a strongly 
renormalized Luttinger liquid behavior for generic filling factors, which should exhibit power-law 
suppression of the conductivity at low temperatures Ef /{kpa) <C ksT <^ Ef, where a is the lattice 
spacing and fcj^a S> 1. Some fully insulating ground-states are expected only for a discrete set of 
integer filling factors for the electronic system. A detailed discussion of the scattering matrix flow 
and its implication for the low energy band structure is given on the example of a square lattice. 

PACS numbers: To be determined 
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I. INTRODUCTION 

For the past two decades, transport properties of quan- 
tum wires have received a lot of attention P, Q . Besides 
metallic systems, two dimensional electron gases induced 
in GaAs/AlGaAs heterostructures have displayed a rich 
variety of quantum effects such as Aharonov-Bohm resis- 
tance oscillations in a ring geometry 0, Q , and persistent 
currents ^ @ . Since the electronic transport mean free 
path in such artificial nanostructures can be as large as 
several micrometers, most scattering processes for elec- 
tronic quasiparticles occur at the nodes between several 
conducting wires. Their influence has been extensively 
studied in the context of the breakdown of the Quantum 
Hall Effect in narrow channels. Experiments have re- 
vealed that the Hall resistance measured in a four probe 
geometry disappears at low magnetic fields 7, 8J. Theo- 
retical studies have emphasized the role of quantum me- 
chanical resonances in the scattering amplitudes of elec- 
trons at the junctions between the main channel and volt- 
age probes llCl|| . In experimental systems, confining 
potentials remain smooth in the vicinity of such junc- 
tions, and this induces a rather robust coUimation mech- 
anism for incoming electrons [ill IT^ . This semi-classical 
picture has been confirmed by spectacular experiments 
involving junctions with various shapes |l3j |. More re- 
cently, coherent Aharonov-Bohm oscillations have been 
measured in ballistic arrays with the dice lattice geom- 
etry [T3 |. in agreement with the predictions of simple 
models for non- interacting electrons [isl IT^ . 

In parallel to this mostly single electron physics, 
dramatic electron-electron interaction effects have been 
demonstrated in transport measurements on various bal- 
listic conductors with very few transverse conduction 
channels. For instance, tunneling into edges of a two- 
dimensional electronic droplet in the Fractional Quan- 



tum Hall Effect (FQHE) regime has shown current ver- 
sus voltage curves with power law behavior (17j in qual- 
itative (though not really quantitative) agreement with 
theoretical models based on the chiral Luttinger liquid 
picture Measurements of shot noise associated to 

tunneling processes from one edge to another through 
a quantum point contact have also provided a convinc- 
ing demonstration of the presence of fractionally charged 
quasiparticles in the FQHE phase 19, 20, 21]. Another 
family of one-dimensional quantum conductors are car- 
bon nanotubes. In particular, single wall nanotubes have 
shown a strong reduction of the single particle density of 
states at low energies 122 ,. 23J . compatible with the Lut- 
tinger liquid model 
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These two main lines of research just outlined can be 
naturally combined and lead us to consider the subject of 
networks of interconnected quantum wires, each of them 
being described as a Luttinger liquid. As a first step in 
this direction, several systems with nanotube crossings 
have been synthesized j2Q, |22|, |2^ 2M ■ On the theoret- 
ical side, many studies of Lut ting er liquids crossing at 
one node are now available |30l l3ll l32l l33l l34l , including 
extensions to more complex geometries |35| . In this pa- 
per, we consider a regular network of Luttinger liquids. 
As already mentioned, the main source of electron scat- 
tering in ballistic structures arises from the nodes of the 
network. For non-interacting electrons, these nodes are 
simply described by a scattering matrix [s^ l37l l3^ , and 
the full band structure (in the absence of disorder) can 
be retrieved from the knowledge of this matrix. However, 
as first shown by Kane and Fisher for a single impurity 
in a Luttinger liquid, interaction effects induce a varia- 
tion of the dressed scattering matrix as a function of the 
incoming electron energy 39, 40J- One way to interpret 
this in physical terms is via the notion of Anderson's or- 
thogonality catastrophe: in the limit of a tunnel barrier. 
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an electron jumping across the barrier leaves a dipolar 
charged excitation which is very far from any eigenstate 
of the interacting system. A rather complicated collec- 
tive relaxation process follows any single electron tunnel- 
ing event. A remarkable prediction made in these works 
is a dramatic qualitative difference between repulsive and 
attractive interactions. In the former case, the effective 
impurity potential grows as the typical energy becomes 
closer to the chemical potential. So a single impurity is 
sufhcient to disconnect completely an infinite Luttinger 
liquid at T = for repulsive interactions. Conversely, 
any static impurity becomes transparent in the low en- 
ergy limit in the case of attractive interactions. 

An appealing picture for these effects has been pro- 
posed by Yue, Glazman and Matveev They have 
shown that renormalization of the transmission ampli- 
tude may be attributed to scattering of an incoming elec- 
tron on Friedel density oscillations induced by the scat- 
terer. From this picture, they have developed an alterna- 
tive renormalization approach, which is perturbative in 
the electron-electron interaction, but non-perturbative in 
the strength of the impurity potential. This framework 
has been used later in references [s^, , and we shall 
adopt a similar procedure here. Note that a third type 
of renormalization scheme, involving the full momentum 
dependence of the electronic self- energy, has been imple- 
mented in a series of papers ■ 

The main novel feature in regular arrays in compari- 
son to simpler geometries as a few connected wires is the 
existence of commensurability effects between the Fermi 
wave-length of electrons and the lattice period. In a non- 
interacting electron picture, we expect an energy gap in 
the spectrum when the average electron number in each 
unit cell of the lattice is an integer. As we shall see 
later, the band structure for a two-dimensional network 
yields a gapped excitation spectrum whenever some in- 
teger numbers of bands are filled. For interacting elec- 
trons, commensurability effects may also be understood 
by considering the pattern of Friedel density oscillations. 
In a one-dimensional geometry, these oscillations exhibit 
a dominant wave- vector equal to 2k p, where kp is the 
Fermi wave-vector for a non-interacting one-dimensional 
wire with the same electronic density. Let us denote by a 
the distance between two nodes. Friedel oscillations origi- 
nating from different nodes share the same global phase if 
2kpa is an integer, which simply means that the average 
number of electrons along a segment of length a is inte- 
ger. Therefore, in the case of repulsive interactions, we 
expect an insulating ground-state in the commensurate 
case, where the Kane-Fisher mechanism will disconnect 
all the wires incoming at the same node. For incommen- 
surate fillings, we predict a strongly renormalized Fermi 
liquid, where the partially filled band crossing the Fermi 
level becomes much less dispersive than for the original 
non-interacting band structure. We suggest that these 
effects should be in principle observable in networks of 
ballistic wires where the electronic density could be con- 
trolled by an uniform gate potential. By changing the 



gate voltage, these systems are expected to undergo a 
succession of metal- insulator transitions. The difference 
between an interacting system and a non-interacting one 
will be manifested by power-law dependences for the con- 
ductance as a function of temperature at fixed bias volt- 
age, or as a function of V at fixed T [39L lioj. provided 
both fcsT and eV remain higher than an energy scale A 
which is the renormalized band splitting in the incom- 
mensurate case, or the single particle gap in the com- 
mensurate one. 

This paper is organized as follows. In sectional we 
consider a simpler problem, namely a one-dimensional 
chain of regularly spaced impurities. We set up a renor- 
malization group method for weakly interacting elec- 
trons which is closely related to those developed in 
Refs m EH, but where the periodicity of the system 
is explicitly taken into account. Section IlIII generalizes 
this approach to any lattice composed of links of the same 
length, assumed to be large compared to the Fermi wave- 
length. We show explicitly that the scattering matrices at 
each node of such lattices is renormalized exactly in the 
same way as for a single node connecting semi-infinite 
wires [S^, Ej . This is the central result of the present 
work. As an illustration with possible experimental rel- 
evance, section IIVI considers a two-dimensional square 
lattice of Luttinger liquids. We show that although the 
evolution of the scattering matrix of the nodes as the 
typical energy scale is reduced yields a rather trivial 
low-energy fixed point where all the links become dis- 
connected, some interesting qualitative changes in the 
quasiparticle band structure take place along this renor- 
malization group flow. 



II. ONE DIMENSIONAL WIRE WITH A 
PERIODIC IMPURITY POTENTIAL 

The goal of this section is to adapt the simple renor- 
malization group procedure initiated in Refs. ,32> .41,j to 
the case of a periodic potential. The main idea developed 
in these works is to dress the bare scattering amplitude 
by a correction due to the interaction of an incoming 
electron with the Friedel density oscillation induced by 
the impurity. This approach treats the electron-electron 
interaction as the perturbation. Because the continuous 
spectrum of particle-hole excitations in the metallic wire 
exhibits a finite density of states down to arbitrary low 
energy, the first order correction to the scattering ampli- 
tude diverges as ln(|fc — kp\d), where k is the incoming 
electron's wave- vector, and d is the spacial range of the 
bare impurity potential. This type of infra-red diver- 
gence is very similar to those encountered in the Kondo 
problem, and Yue et al. proposed to treat them with 
a renormalization group method inspired by Anderson's 
"poor man's scaling^^ approach [45|. The idea is to inte- 
grate out gradually single particle-hole excitations which 
participate in the Friedel oscillation, starting from those 
furthest from the Fermi level. As the electron bandwidth 
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D is continuously reduced, the bare impurity potential is 
renormalized so that the low energy physical properties 
of the system are kept unchanged. The renormalization 
procedure stops at a low energy scale with is the larger 
scale among the thermal broadening /csT, the bias volt- 
age eV, or the incoming electron energy h\k — kF\vF- 

As already stated in the Introduction, the presence of 
an array of scattering centers (such as nodes in a wire 
network) brings qualitatively new features. In the low 
energy regime, Friedel oscillations originating from dif- 
ferent centers are expected to interfere, so we cannot fol- 
low the renormalization flow obtained in Refs for 
a single scatterer down to arbitrary low energies. Fur- 
thermore, commensuration effects between the average 
electronic density and the superlattice structure play a 
crucial role. By contrast to the single impurity case, we 
expect an insulating ground-state only for an integer av- 
erage filling of each supercell. For incommensurate filling 
factors, we expect a cross-over from the one-dimensional 
behavior following Kane and Fisher's predictions at high 
energy, towards a strongly renormalized coherent conduc- 
tor at low energy with finite conductance. In a periodic 
system, the natural way to implement this "poor man's 
scalin(f^ approach is to integrate out energy bands one af- 
ter the other, starting from those most remote from the 
Fermi level. In the incommensurate case, the last band, 
which crosses the Fermi level is partially filled, so it is nat- 
ural to stop the procedure after the last fully occupied 
band has been integrated out. In any renormalization 
method, we have to decide which low-energy quantities 
will be required to remain constant as high energy modes 
are eliminated. In the presence of a periodic potential, it 
is natural to prescribe that single quasiparticle energies 
should not change under the renormalization group flow 
(RGF). 




FIG. 1: Localized impurity potential may be represented by a 
5'-matrix, that connects amplitudes of incoming {A, B') and 
outgoing {A',B) plane waves outside the impurity. 



superposition of plane-waves: 



for X < —d 
ioT X > d 



(2) 



A. Band structure for a periodic array of point 
scatterers 

Let us first consider a non-interacting problem along an 
infinite one-dimensional wire with a periodic potential. 
The corresponding Hamiltonian is 



-2 

H = Vix-na) 

n=0 



(1) 



where a denotes the spacial period of the potential, 
namely the distance between two succesive impurities. 
V{x) is a localized potential, so for instance we impose 
that V{x) = when |a;| is larger than a range d, d <^ a. 
The effect of each scatterer is described by a scatter- 
ing matrix S. Suppose first we have only one of them, 
centered at the origin a; = 0. Let us consider scatter- 
ing states with the energy Eo{k) = fi,^fc^/(2m), k being 
positive. Away from the impurity, that is for > d, 
we may represent the corresponding wave-function as a 



Since Schrodinger's equation is linear and of second 
order, we may express the outgoing amplitudes A' and 
B linearly as a function of the incoming ones A and B' 



A' 
B 



t r' 
r t' 



A 
B' 



S 



A 
B' 



(3) 



where {r, t, r', t'} are two pairs of reflection and transmis- 
sion coefhcients for left and right sides of the node. In 
principle, these four coefficients do depend on the energy 
of the particle or equivalently on its wave-vector k. In 
this paper, we shall neglect this variation, since the dom- 
inant contribution processes involve virtual excitation of 
particle-hole pairs in the vicinity of the Fermi level. A 
more complete approach would consider the Taylor ex- 
pansion of S in powers oi k — kp, but all terms beyond 
the O'th order one are irrelevant according to the clas- 
sification of perturbations around a non-interacting one- 
dimensional fermion system. At least for not too large 
interactions, they are not supposed to change the qual- 
itative picture of the system behavior. As usual, this S 
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matrix is unitary. Assuming time reversal invariance of 
the Hamiltonian implies t = t' and if V{x) is an even 
function of x, we have also r = r' . In this case, we may 
parametrize S by two angles {0 < (f> < 7r/2, < '0 < 27r) 

g ^ f cos (j) ±i sin cj) \ 
I ±z sin (j) cos y 

For a periodic array of identical scatterers, eigenstates 
may be obtained as Bloch functions, namely we may im- 
pose the condition 

where k' is chosen in the first Brillouin zone [— 7r/a, tt/o]. 
On each x-interval [na + d, (n + l)a — d], we write the 
eigenstate with energy Ea{k) as 

The above periodicity condition implies 

A — J(k'-k)anA 

Eq. (jSJ can now be written for each impurity site, which 
gives: 

(5) 

Replacing 74„ and B„ by their expressions in terms of 
and Bq, we get the following secular equation 

^gi2fea ^l^i(k+k')a _ ^ ^ (^) 

which determines the dispersion relation implicitly via 
fc, the energy being EQ{k), the lattice momentum k' be- 
having as an external parameter. Using a normalization 
condition on the wave-function, we could get {Aq, Bq) as 
functions of {k' , S). 

In the particular case of spacially even and time- 
reversal invariant potentials, we may use the above 
parametrization for S in Eqn. ( |HJl , which yields 

cos(fca + ip) = cos (j) cos(fc'a) (7) 

For a given value of the lattice momentum fc', the possible 
values of ka appear in two equally spaced families, with a 



period 27r for each of them. The allowed values of ka + ip 
belong to the intervals [— tt + + 27rn, —0 -I- 27m] and 
[(j) + 27m, TT — (j) + 27m], where n is integer. We recall 
that k should be positive, in order not to count each 
eigenstate twice. The values of ka + ^ lying in intervals 
[— (/) + 27rn, + 27rn], and [t: — 4> + 27rn,7r + (f) + 27rn] 
correspond then to energy gaps. These gaps are of course 
larger when the reflexion coefficient is larger. 




FIG. 2: Band structure for ID wire of noninteracting electrons 
with periodic impurities. BS is 2n periodic in k. A set of gaps 
{ — (f> + 2nn,<l) + 2nn) and (tt — + 2?™, tt + + 27rn) is present 
for any value of S-matrix. 



B. Switching on electron-electron interactions 



In a Luttinger liquid, the effective interaction becomes 
non-local in the low-energy limit. To show this, it is 
convenient to decompose the electron creation operators 
(where a G {T,J.} denotes the spin component 
along a fixed direction) into a right moving part '^~j^„{x) 
and a left-moving part 4'Jcr(^)i where "^^^ (resp. ^'J^.) 
involves the Fourier modes k close to kp (resp. —kp)- 
With this decomposition, the local electron density p{x) 
is written as follows 



p{x)= vl/+(a:)vl/,(x) - + E (*L(^)*i^'^W + *LW*iia(x)) (8) 

The first two terms are smooth fields, meaning that their Fourier transforms involve only small wave- vectors compared 
to kp. But the last two terms are centered around the wave-vectors ±.2kp so they are rapidly oscillating. 
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For a spin-rotation invariant Hamiltonian, the effective 
low energy description of a Luttinger liquid involves three 
independent parameters: the velocities Vc and Vg of col- 
lective charge and spin excitations, and a dimensionless 
constant K which depends on the strength of electron- 
electron interactions and controls the exponents enter- 
ing the correlation functions. Since transport proper- 
ties are mostly affected by the value of K [33, we shall 
not consider here the renormalizations of Vc and Vg away 
from their common value vp for a non-interacting sys- 
tem. Therefore, it is sufhcient to consider the following 
interaction 



Hint = —r- 



2 



L/2 



(9) 



-L/2 



where po (x) is the long wave-length part of the total den- 
sity: 

Po{x)= ^ {^+^{x)^nA^) + HA^)^LA^)) 

Here L denotes the total length of the system. Later, we 
shall assume periodic boundary conditions, and that L 
encloses an integer number N of periodic cells, so L = 
Na. With this choice of interaction, we have 



Vc = VF 
Vs = Vp 

K ^ (1 



1/2 



2Uo_\ 



-1/2 



So ii" = 1 for a non-interacting system, K > 1 for attrac- 
tive interactions, and K < 1 for repulsive interactions. 
For our purpose, it is convenient to view this effective in- 
teraction as deriving from a non-local potential U{x — y) 
such that its Fourier transform U{k) vanishes outside a 
finite window centered around fc = and whose width is 
smaller than 2kp. The interaction strength Uq is defined 
as U{k = 0). With this notation, we have 



L/2 L/2 

J7int = i j dx j dypQ{x)U{x 

-L/2 -L/2 



y)po{y) (10) 



In this section, we are considering the combined ef- 
fect of impurity scattering and interactions. Renormal- 
izations of the effective scattering matrix S are natu- 
rally detected via the electron self-energy S(fc, k' , uj). But 
since our system exhibits only a discrete translation sym- 
metry, we may only conclude that k' — k should be an 
integer multiple of the basic reciprocal lattice vector 
This self-energy is then a relatively complicated object. 
More information on its rea l-sp ace structure for a single 
impurity may be found in |43 . l43j |. To analyze it in a 
simple way we shall compute the first order correction 
Ei{k) with respect to Uq to the single electron energy 
Eo{k) = h^k'^ / {2m). Here k stands for a single parti- 
cle level, close to the Fermi energy, and labelled by the 
combination of a Bloch quasi-momentum k' and a band 
index. This correction Eiik) is given by the sum of a 
Hartree term and of an exchange term. In the case of an 
unpolarised electron system, we have 



L/2 L/2 

Ei{k)= f dx f dyJ2 rk{xWg{y)U{x-y)[2i;g{y)Mx)-M^)My)] 



(11) 



-L/2 -L/2 



q<kf 



For local potentials, the Hartree and the exchange con- 
tributions cancel each other, when the spins of the two 
electrons involved are parallel. But as we have recalled 
before, our two-body effective potential is in fact non- 
local, so we have to analyze both terms in more detail. 
Our expression for Ei{k) involves integrals of the form: 

L/2 L/2 

Iu{f,g)^ J dx J dyf*{x)U{x ~ y)g{y) 

-L/2 -L/2 

where f{x) and g{y) are Bloch functions satisfying 
— ff X = r-— = e , < 6^ < 27r 



Since U is short-ranged in space (though it is not a delta 
function), we may take safely the thermodynamic limit 
i ^ 00. Writing f{x) = /,^e^«2^"+^)^/'^) and an 
analogous series for g{y), we obtain: 



/^(/,g) = L^/,:C/(^^^)5„^LC/o ^ /^ffn 

» |n|<fcj?a 

(12) 

Let us first consider the Hartree term. Because we 
chose single particle eigenstates of the Bloch form, the 
corresponding local particle density is periodic with pe- 
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riod a. A little elementary algebra shows that 

2 1 sin(fca + '0) ± sin(0) cos(2fc(a; — f )) 

h/'fci'^) ~ — 

Na sin(A:a + t/)) ± sin((/)) sin(/ca)/(fca) 

(13) 

for < a; < a. As expected, the amplitude of the local 
density oscillation is stronger when the bare reflexion co- 
efficient is larger, or equivalently when | sin(0)| is larger. 
The n}^ Fourier amplitude of this local density is equal 
(for n^Q) to: 

Ak ( sin(A:a — 7rn) sin(fca + tto) \ 
IN a \ ka — irn ka + nn J 

where the numerical coefficient Ak is close to unity. As 
shown in Eq. H12|l above, we are interested in the case 



where |n| <C fc^a, and since k is close to fcf, this ampli- 
tude is small by a factor l/(kpa). A similar conclusion 
holds for the Fourier amplitudes of |4'q(a;)p if we assume 
that the most important effects come from filled states 
where q is close to fci?- Therefore, we do not expect strong 
renormalizations coming from the Hartree term. 



Let us now turn to the exchange term. The product 
^*q{y)^k{y) is the sum of four oscillating terms propor- 
tional to e^'C^^^)^ and e±*('=+«)''. The last two terms are 
fast oscillations which will be filtered out by the non-local 
potential, as in Ea. H12|l . Keeping only the first two oscil- 
lations, we can cast the exchange contribution to Ei{k) 
as follows: 



/, \ t^o /,/ ,T N Uq . y ^ f dq' sm(q(q') ~ k) , ^, 

El (k) = c— (fcp + Nftt) shr 6 / ^ , , , — ^ , ^ (14) 

^ ' a ^ ^ ' 27ra siu{k + ip) shi{q{q') + ^j) q{q') - k ^ ' 

In this equation, we have replaced combinations such as ka, qa, by new dimensionless variables fc, q. The integral 
symbol stands for a summation over all the Np completely filled bands, including possibly a last partially filled band 
with a dimensionless momentum k'p such that < fc^ < tt. For each completely filled band, the integration variable 
q' runs from to tt, and q in Eq. H14() is a function of the lattice momentum q' solution of the dispersion relation ((T)). 
For the last partially filled band (incommensurate case), the q' integral runs from to k'p. As already mentioned, we 
have assumed that parameters {ip, (j)) are not depending on the incoming energy. Note that contributions from the 
Hartree term will modify only the numerical coefficient c whose precise value is not important here. 



C. Renormalization approach 

Let us introduce the notation Ag = ttNp, which plays the role of a large momentum cut-off. As in all schemes 
inspired by Anderson's "poor man's scaling", we shall assume it is possible to construct a sequence of models where 
filled bands are eliminated one after the other, starting from the most remote from the Fermi level. When the first n 
bands have been eliminated, the new value of A is set equal to Aq — nn. At each step, we require that the quasiparticle 
energy [Etot{k) = Eo{k) + Ei{k) = Ea{k) -f Uo£i{k)) for k close to the Fermi wave- vector should remain unchanged. 
To compensate for the reduction of the cut-off from Aq to A, we have to adjust S'-matrix parameters {"0, </>} so they 
become functions of running cutoff A. This is expressed by the following prescription 

Etoti-^o. 00, k) = So(^(A), 0(A), fc) + C/o£i(V^(A), 0(A), fc, A) (15) 

Since for Uq = 0, this condition implies i5o(V'Oj 0o, fc) = £'o(V'(A), 0(A), fc), we see that in this case (?/;o,0o) = 
('0(A), 0(A)) for any A, so we may write the following Taylor series 

^(A,C/o,^/>o,0o,Ao) = Vo + C/o?(A,Vo,0o,Ao)-l-O([/o^) 
0(A,f/o,^/'o,0o,Ao) = 0o + C/o0(A,0o,0o,Ao) + O(C/o) 

We now try to keep band structure (|15(l unchanged for any k 

Eoii^o, 00, A:) + UoSiii^o, 0o, fc, Aq) = £;o(V'o + W(A), 0o + C/o0(A), /c) + C/ofi(?Ao + C/oV'(A), 0o + C/o0(A), /c, A) 

Keeping the first order terms in Uo gives 



dEo 



dip 



n 771 

^(^) + ^ 0(A) = 5i(V'o,0o,fc,Ao)-^i(0^o,0o,fc,A) (16) 

V-o O0 ' 



r 



This is a non trivial constraint, since ^^\'4>o a-nd ^^l^o "0 and depend on A but not on k. The possibility to 
depend on k but do not depend on A. On the contrary. 
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enforce this requirement is not obvious a priori, and when 
it occurs, we may call our model renormalizahle (at least 
to this lowest order). 

Let us now evaluate the right-hand side of this equa- 
tion. Suppose we integrate out just one band, then 
Aq — A = vr, which is assumed to be much smaller than 
A. While computing £i(V'o, 0Oi Aq) — £i{ipQ,4)Q,k, K) 
in Eq. (|14|l . the integral involves only one band far from 
the Fermi level. Therefore, we may further approximate 
q{q') - fc by -A. This yields 



£:i(V'o,0o,fc,Ao) - £i(V'o, 00, fc, A) 



c(Ao - A) 



1 1 



sm 



r , , sin(g(gO - fc) 



A 27ra sin(fc -I- ^o) Jo sin(g(g') -f t/jq) 

From the derivatives involved in the left hand-side of 
Eq. l(T^ are: 



dtp 
dEo 



(p.k— const a 

« tan(pcot(K -I- ?/') 

i/i-A:— const ft 



We have linearized the bare dispersion relation: Eo{k) 



2rn(P 



const ■ 



k. The notation "fc = const" means 



more precisely that the Bloch crystal momentum k' and 
the band index have to be maintained constant while 
varying (j) or tp. Introducing these expressions for the 
derivatives and the result H17|l in Eq. H16|l shows that 
indeed the k dependences on both sides can be made 
to match, which expresses the renormalizability of our 
model to first order in interaction strength. This fixes 
the form of the functions V'(A) and (/'(A): 



V'(A,'(/'o,Ao) = 



</>(A,0o,Ao) 



!^lnA/,,'eot(,(,') 



^ii; ^^"(200) In A 



(po + V'o)) 



Finally we construct the RGF equation: 



^° fc-H^i /d,'cot(,(,')+^)) 



dA nhvp 27r A 







Un 



9 In A A-Khvp 



sin 20 



(18) 



We see from Eq. (0J that the parameter -0 is a global 
phase in the scattering matrix, which does not affect any 
physical property of the system besides an overall shift 
of the single particle spectrum. In particular, it does not 
generate any density oscillation. Moreover the associated 
RGF equation explicitly involves the running cut-off A, 
and the notion of fixed point loses its meaning here. 



Therefore, we now turn to 0(A), for which a simple 
RGF equation arises, and which solution is given by: 



tan0= (Ao/A)"tan0o 



(19) 



where a = Uq/ (2TThvp). The corresponding transmission 
coefficient r(A) on a given impurity is: 



T(A) = cos2(,/)(A)) = 



ro(A/Ao)' 



i?o+To(A/Ao)2 



(20) 



where Tq is the transmission coefficient for a single im- 
purity in absence of interaction, and Rq — 1 — Tq. This 
result agrees with the expression obtained for a single im- 
purity 41] in the absence of spin backscattering, namely 
when U{2kp) = 0. Again, this approach assumes small 
electron-electron interactions. In the case of strong in- 
teractions, where K is no longer close to 1, the bosoniza- 
tion method shows that for the single impurity problem, 
a should be replaced by (1 — K)/2 These two ex- 
pressions for the exponent coincide at small Uq if terms 
of order Uq or higher are neglected. 

For a commensurate system {kpa — nn, n integer), the 
non-interacting ground-state is already gapped, so we ex- 
pect a true insulator as well in the presence of repulsive 
interactions. The difference between a traditional band 
insulator and the one obtained here in the presence of in- 
teractions is the non-trivial energy dependence of the im- 
purity scattering matrix and the corresponding behavior 
of the Landauer conductance, Eq. (|20|l . For an incom- 
mensurate system, we have a partially filled band cross- 
ing the Fermi level in the absence of interaction. Since 
our renormalization procedure assumed a gradual elimi- 
nation of fully occupied bands, it has to break down after 
the last of those bands has been integrated out. Treat- 
ing the remaining partially occupied band in a heuristic 
way, we simply assume that it corresponds to a strongly 
renormalized Luttinger liquid, whose effective Fermi ve- 
locity Vp is much reduced compared to the Fermi velocity 
vp — hkp/m of an uniform non-interacting gas with the 
same density. More precisely, we have: 

Vp ~ Vp cos 05 sin(fc^a) 

where 0s ~ 7r/2 is the value of when the renormaliza- 
tion procedure stops, which corresponds to: 

Aq kpa 

A TT 



Using Eq. (fll^ . we get: 



Vp = Vp \ I cot 00 sin(fcpa) 

kpa : 



(21) 



III. 



GENERALIZATION OF RG PROCEDURE 
TO A LARGE CLASS OF LATTICES 



In the previous section we introduced the main ideas 
we used to obtain the RGF equation for a ID lattice. We 
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wish now to show that renormahzabihty of this particu- 
lar ID system is not a simple coincidence, but a general 
property of any network (not necessarily periodic), pro- 
vided the two following assumptions hold, namely all the 
links have the same length, which has to be large com- 
pared to the Fermi wave-length. Let us begin to follow 
the same procedure as in one dimension. Suppose that we 
have a network of equal length wires. Any junction point 
is described by an unitary ^-matrix which dimension is 
equal to the number of wires joining at this node. For 
each link, stationary single electron states can be written 
as the sum of two plane waves. 

^Pix) = A.jC-'''^ + Aj.e''''' (22) 

where Aij is the amplitude of the wave that propagates 
from node j to node i, if the x coordinate is oriented 
from i to j. Solving Schrodinger's equation is equivalent 
to connect these various amplitudes via node scattering 
matrices: 

A,^Y.^'\'"'si^2A,m (23) 

m 

Here means that we sum over first neighbors m of 

m 

node j. We notice that this has indeed the form of an 
eigenvalue equation written in some basis. Following the 
idea of Kottos and Smilansky we introduce a finite 
dimensional Hilbert space associated to the lattice links. 
Each link ij is represented by two orthonormal vectors 
and \ji). The dimension of this auxiliary Hilbert 
space is therefore 2Nl {Nl is the total number of links). 
One may rewrite Ea. H23|) in its vector form: 

f\A{k))=e-'''''\A{k)) (24) 

where the T operator incorporates information about the 
scattering matrices of all nodes. 

^ = E E''^ mS^/2{jm\ 4^ = {^J\f\Jm,) (25) 

j i,m 

I 



As this operator T is unitary and defined in a finite 
dimensional Hilbert space, it could be diagonalized as: 
T\a) — e^'^°|a), where a takes 2Nl values and 9a is 
real. So we obtain families of eigenvalues for the single 
electron energy E — h'^k'^ /2m: 

aka,„ ^0a + 2Tm > (26) 

We emphasize that this periodic structure of the single 
particle spectrum is a special feature of constant link 
length networks. A brief discussion of the more general 
case is given in appendix El Because of this periodicity, 
and despite the absence of any translational symmetry, 
we may still introduce a notion of energy band for such 
lattices. More precisely, in this setting, an energy band 
corresponds to fixing n and allowing for all possible val- 
ues of 9. Note that this notion of band does not exactly 
coincide with the more familiar notion from the Bloch 
theory of translational invariant lattices. For simple Bra- 
vais lattices, the number of states in each Bloch band is 
the number of unit cells which is equal to the number 
of sites Ns- If Z is the coordination number, we have 
ZNs = 2Nl, so our generalized bands contain Z usual 
Bloch bands for a Bravais lattice. At this stage, we have 
so far a band structure equation written in operator form. 
In order to obtain renormalization flow for 5-matrix we 
need to compute first the electron-electron contribution 
as in Eq. Hll|) to the single electron energy and then the 
variation in the unperturbed energy due to an arbitrary 
^-matrix change dE^/dS. 

As in one dimension, the main contribution to the elec- 
tronic self-energy is given by the exchange term. Let us 
consider a pair of single particle eigenstates labelled by 
k and g, where these labels should in fact be viewed as 
pairs (a, n) and (/3, m), m and n being integers according 
to the above description of the spectrum. State k is close 
to the Fermi level, but state q is far from it, at a distance 
corresponding to the current energy cut-off A. Along a 
link ij, we denote by (\l/^(a::)\E'q(a;))o the slowly varying 
component of '^^.{x)^' q{x) . A simple computation shows 
that: 



|(vI/*(x)*,(x))o|^ = \A,mM^A<l)? + \A,mM,^{<l)? + (^*.(fc)A,,(fc)A,,(g)A^(<z) +h.c.)^^|^^-^(27) 
The first part summed over fully completed band does not depend on energy: 

Y,\A^,{k)\M^A<l)? + \Arm?\Ad<l)\''= E E = 

g^Band i,j q^Band <ij> 

= {i]\k){k\ij){i]\ Y - E = = 1 (28) 

<ij> qgBand <ij> 

We used the fact that both \q) and \ij) form complete basis sets in our Hilbert space. The main expression to compute 
is then 

E E ^j™(^)^"^(^)^^™(9)^™j(9)sin(/c-g)a (29) 

(jGBand <jm> 
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where the sum over q is just a single band sum, namely 
m is fixed, and the sum is taken over the 2N]^ values 
of (3. The power of this algebraic formalism is that such 
sum is readily performed, without having to compute any 
integral. Indeed, we have: 

^(^) E {m){^MM[k^q)a]{q\ji){ji\k) (30) 

q <ij> 

As shown in appendix [HI we may assume that the 
eigenvectors \q) are normalized to unity in the auxiliary 
Hilbert space attached to link amplitudes, provided the 
links have the same length, much larger than the Fermi 
wave-length. Therefore, we have the very useful com- 
pleteness relation, that is: 



(31) 



After some simple algebra, we may cast I{k) into the 
form 



/(fc) 



g»fca_^^|Y|^^ (32) 



where the single node operators V^'^-* are defined by: 

yU) = pU) _ gij) p{3)tgU) 

and the diagonal matrix i^*^^^ by: 



We have introduced as before the dimensionless parame- 
ter a = Uo/{2Trhvp). 

To get the first order variation of the single electron 
energy under small changes in the node scattering matrix 
parameters we differentiate Eq. (|24|) : 

idf)\k) + f\dk) = ~i{dk)ae-''^''\k) e-'^''\dk) (33) 

Applying (fc| to this equation and using {k\T = e^''^°(fc| 
we obtain: 



{k\k) 



(34) 



This allows us to calculate single electron energy varia- 
tions due to 5'-matrix changes. In the particular case 
of global phase transformation, the corresponding in- 
finitesimal form reads: dT = iTdip. Clearly, the en- 
ergy differential does not depend on energy any more 
since dk = —dtp /a, so global phase shifts simply induce 
a global translation on the energy spectrum. 

Following the same ideas as in ID, we generalize the 
RGF equation to any S'-matrix parametrization. Eq. H16|l 
now becomes: 



dEp 
dS 



{So,k)dS{A) ^ E,{So,k,A„) - E,{S„,k,A) (35) 

The left-hand side of this equation is equal to hvpdk, 
where dk is related to the small renormalization of T 
by Eq. To evaluate right-hand side, as before, we 

integrate just over one band of width 2tt for the quantity 
Qp.rnO-, i-e. dA = A — Ao = — 27r. 



nvpdk = Ei{So,k,Ao) ~ Ei{So,k,A) = 

A<qa<Ao <jm> \ 1/ \ / 



Constant c includes both the Hartree term and the part 
of exchange term that does not depend on k, so we do 
not precise its value since it renormalizes only the global 
phase of S'-matrix. 

Finally, we get the result dT/dl = —V that agrees com- 
pletely with Lai, Rao and Sen [^3, obtained for a single 
node connecting an arbitrary number of scmi-infinitc ID 
wires. In coordinate way of writing, it gives: 

^ = s(^")F(^")ts(^) - f(^) (36) 

where we just chose the usual cutoff parametrization: 
A = Aoe-'. 



IV. TWO-DIMENSIONAL SQUARE LATTICE 



We would like to illustrate the result of the previous 
section on one more example. This part could be in- 
teresting from an experimental viewpoint, since present 
nanofabrication techniques are now available to prepare 
networks of quantum wires with a very small num- 
ber of transverse conduction channels etched on a two- 
dimensional electron gas with high mobility, as illustrated 
for instance in |14|. Let us now consider an infinite reg- 
ular square lattice of perfect Luttinger wires. These one 
dimensional conductors are only coupled at the lattice 
nodes which are described by a single 4x4 scattering ma- 
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ft., 



FIG. 3: Two dimensional periodic grid of electron liquids with 
impurities. Each impurity could be represented by 3 complex 
parameters: r- reflection, ty- forward transmission, t±- per- 
pendicular transmission coefficients. 



trix S. To keep a simple model, we shall restrict ourselves 
to the case of a single conduction channel in each wire, 
although the case of several channels would clearly be of 
interest, both on the theoretical side, and with respect 
to possible experimental realizations. As mentioned in 
the Introduction, we shall not take into account any en- 
ergy dependence of the scattering matrix, although de- 
tailed studies of the Schrodinger equation for a cross of 
wires with a finite width have exhibited a rich pattern of 
resonances [4^ l48j . The main motivation for this simpli- 
fied treatment is that in a renormalization group picture, 
smooth energy dependences in the scattering matrix as 
a function of E — Ep correspond to irrelevant operators, 
which should not alter drastically the way interactions 
drive the system to its low-energy fixed point. Labelling 
the four directions joining at a node as on Fig.|21 we shall 
consider a scattering matrix of the following form: 




(37) 



which corresponds to the most general form obeying time 
inversion and spacial dihedral symmetry, in combina- 
tion with unitarity. The previous expressions involves 
three complex parameters, but as shown in appendix 1X1 
unitarity leaves only three independent real variables. 
We have chosen the following parametrization: 



2)/4 

t|l = e^'l' {e^i't'- +e^"t'" +2) /A 



(38) 



)/4 



where (t)u,v G [0, 7r[ and V' £ [0, 2tt[ 
Note that two lines in the {(pu,'^ 
interesting: 

f (f>u = 7r/2 ^ t_i_= t|| 

1 (I)V. = tl = 



plane are specially 



(symmetric case) 
(ID case) 



A. Band Structure 

The derivation of the band structure is standard, so it 
is outlined in appendix IbI This band structure is given 
by an implicit equation: 



x(fc,k')+2/(fc,k') =/3 



2 cos ( 



sin(^„ - (j)y) 



where 
x(fc,k') 



sin(fca + tp) 



cos(f>u cos(fc^a) — cos(fca + ip + </'«) 
,^ _ sin(fca + ip) 

cos(j)u cos(fc' a) — cos(fca + ^p + 4>u) 



(40) 

(41) 
(42) 



As usual, the energy of these states is given by the 
free electron dispersion Eo{k) — ft^fc^/(2ma^). Here, 
k' is the two-dimensional lattice wave-vector, such that 
^'(r-l-R) — exp(ik'.R)^'(r) for any r on the wire net- 
work and any period R of the square lattice. 

As we found some interesting features in the band 
structure of noninteracting electrons in a two dimensional 
square grid we will describe it more precisely. Contrary 
to one dimension, there are values of the scattering ma- 
trix, for which the single electron spectrum is no longer 
gapped, and these are located on Fig. 0] More precisely. 




(39) 



FIG. 4: Phase diagram for a non-interacting electron wire 
square grid. Contrary to the one dimensional case, there are 
metallic states at integer filling factors for some values of the 
S-matrix. In these regions of the phase diagram, the single 
electron spectrum is gap less. 

in the clear regions of Fig.^ the single particle spectrum 
is gapless. In the dark regions, it is gapped, leading to 
an insulator if the electronic density corresponds to filling 
an even integer number of bands. Finally, in the dashed 
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regions, we obtain an insulator for an odd integer number 
of bands. 

We still have a 2t: periodic structure in fca, but the 
band-structure consists of two types of foils: normal and 
abnormal. Normal bands resemble an ordinary band of 
a tight-binding model of square lattice crystal (a sort of 
deformed paraboloid). Abnormal bands are so called for 
their strange curvature. To get an idea of their form one 
could imagine a square rubber foil, attach its four extrem- 
ities and then put inside a heavy cross. For a complete 
description, we give sections of the band structure in sev- 
eral directions for three characteristic values of the scat- 
tering matrix. Because of some important symmetries, 
we may restrict the domain of variation of {</)«, 0d}, and 
still get all the possible different physical pictures: 



1. k{4>u,4'v) 

2. fc(7r — 01,, 



k{4'v,4>u), 

- 0«) = -k{4>u, 4>v)- 



These may be easily seen from form II of the dispersion 
relation, given in appendix B. Both of them are reflection 
symmetries. Given the band-structure for k G [0, tt] and 
using the following symmetry: k{k'^,k'y) + n = kik'^ + 
TT, k'y -l- tt), we easily expand it to the full interval k G 
[0, 27r] by replotting the same band originating from point 
M instead of T (see Fig. E)). 



B. RGF equation for a two dimensional grid 

Following the same procedure as in the one dimen- 
sional case, we first calculate Hartree and exchange con- 
tributions to single electron energy and then establish the 
equivalent of Eq. H16|) or Eq. H35|l for two dimensions and 
finally get the RGF equation. The main difference with 
the ID case is that we have now three real parameters for 
the S-matrix and electron-electron interactions should be 
evaluated along two perpendicular threads that form our 
grid. The condition to satisfy now reads: 



dEo 



dEo- 



= £i(5o,k',Ao)-£i(5o,k',A) 



(43) 



As was proven in the previous section all networks with 
links of equal length are renormalizable i.e. there is a 
set of functions 4'ut4'v ^^'^ ^ depending only on A. In- 
deed the decomposition of the r.h.s. of Eq. H43|l on a 
basis of three functions depending on k' is possible. The 
corresponding renormalization group flow equations are: 



d<j)u 
dl 



-(sin2(?!)„ + 3sin20„ + sin2(0„ - 0„)) (44) 



^ = ^(sin2(?!)„ -f 3sin20,; -I- sin2(0„ - 01,)) (45) 

Ml O 

where a is defined in Eq. (|40I41II . The only fixed points 
are (j)u,v = Oi 7i'/2, among which there is only one attrac- 
tor for {0u,0t,} = {7r/2,7r/2}. The global behavior of 




b) 
























-C) 

























FIG. 5; 3 characteristic band structure pictures for different 
value of S'-matrix: 

(a) Insulator, < (f>u,<i>v < 7r/2, (dark on Fig|lJ 

(b) Insulator, <<?!)„ < 7r/2,7r/2 < 0^ < vr, |0„ - 0„j < tv/2, 
(dashed on Fig^J 

(c) Conductor, < 4>u < 7r/2, tv/2 < 4>^ < tt, \4>u — 4>v\ > tv/2, 
(clear on Fig|lJ 

The band-structure is 2tv periodic in k, and has 4 foils:2 nor- 
mal and 2 abnormal. Some foils are described as "abnormal" 
because of their strange curvature, revealed here by the flat 
part of these bands. Given the energy interval < fc < vr one 
could obtain the tv < k <2tv interval by exchanging F and M 
points. 



this flow is illustrated on Fig. These properties of the 
RGF for a single node connecting four semi-infinite wires 
have already been described by Lai et al. [s^ and Das 
et al. 35]. As for the one dimensional example of sec- 
tion^ above, the new feature associated to a regular lat- 
tice is the presence of commensurability effects. We have 
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FIG. 6: RGF for 2D square grid of wire. The only attractor 
is in the picture center for {4>u, (pv} = ('""/S, 7r/2). 



to stop the renormalization procedure when all the com- 
pletely filled bands have been eliminated. From Fig. [3 we 
expect to obtain one or two partially filled bands crossing 
the Fermi level. These bands are only very weakly dis- 
persive, since the effective S matrix for the nodes is then 
very close to its value at the vanishing transmission fixed 
point. Suppose now that this fixed point is approached 
from the dark regions of the phase-diagram shown on 
Fig. 01 If the filling factor corresponds to an even in- 
teger, the Fermi level lies in a gap of the renormalized 
band structure. Therefore, we may eliminate the remain- 
ing pair of filled bands, and the system is an insulator. 
Similarly, a true insulator is obtained for an odd integer 
filling factor, in the case where the (7r/2, 7r/2) fixed point 
is approached from the dashed regions in Fig. 0] Ex- 
perimentally, one expects transitions between these com- 
mensurate insulators and strongly renormalized "heavy 
electron" metals at generic filling factors if the electronic 
density is controlled by a uniform external gate voltage. 

Another interesting feature of this geometry is the fact 
that the fiow may induce metal-insulator transitions for 
some commensurate filling factors at a finite energy scale. 
Indeed, for initial parameters lying in the clear regions 
of Fig. 21 corresponding to a gapless single electron spec- 
trum. Fig. (HI shows that the system always reaches either 
the dashed or dark regions in a finite RG time. Experi- 
mentally, these RG fiows may be visualized by gradually 
lowering the temperature, since at least qualitatively, the 
energy scale set by temperature plays the role of the mov- 
ing cut-off A. 



V. CONCLUSION 

In this paper, we have studied a particular class of 
networks of Luttinger liquids, with nodes connected by 
links of a constant length. In the limit of long links, com- 
pared to the Fermi wave-length, we studied the evolution 
of the scattering matrix at the nodes, as the typical en- 
ergy scale for the occupied states contributing to Friedel 
oscillations is getting closer to the Fermi level. The cor- 
responding renormalization group fiow turns to be iden- 
tical to the one already found for a single node coupled 
to several semi-infinite ID Luttinger liquids [s^. This 
result is physically reasonable, since we have considered 
the limit of long links. However, we emphasize that these 
renormalization effects come from quasiparticle scatter- 
ing on Friedel oscillations induced by the nodes, which 
are a rather complicated function of the lattice geometry. 
For instance, even in the limit of very long links, the am- 
plitudes Aij which determine the value of energy eigen- 
functions along the links are obtained from a 2Nl x 2Nl 
eigenvalue problem whose solution has a strongly non- 
local character. 

The main difference between a regular lattice and a 
simple node coupled to infinite wires is that in the for- 
mer case, we have to stop the renormalization procedure 
when the last occupied band has been integrated out. So 
instead of having completely disconnected wires in the 
low-energy limit, we expect in general a strongly renor- 
malized conducting system with an effective Fermi veloc- 
ity much reduced in comparison to a non-interacting sys- 
tem with the same density. These effects should be visi- 
ble as a power-law behavior of the network conductance 
as a function of temperature. Insulating ground-states 
are expected when the electronic density corresponds to 
filling some integer numbers of bands. 

Of course, this work leaves many open questions. It 
would be interesting to generalize the present renormal- 
ization approach to lattices containing links with several 
different lengths. In such situations, the spectrum no 
longer exhibits a simple periodic structure, and some sig- 
natures of quantum chaos, already manifested in the sin- 
gle particle density of states , may also appear in the 
temperature dependence of the conductivity of an inter- 
acting system. Another open question is the influence of 
an external magnetic field, which also drastically modifies 
the single-particle spectrum. Finally, the limit of strong 
electron-electron interaction deserves further investiga- 
tion, and in particular the possibility to develop some 
new metal-insulator transitions for non-integer but ra- 
tional filling factors, generalizing the notion of a Wigner 
crystal. Such insulating states would naturally be pinned 
by the nodes of the lattice. 
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on various aspects of conducting wire networks. 



APPENDIX A: PARAMETERIZATION OF THE S 
MATRIX 

In this appendix we will show that time-inversion, spa- 
cial D4 dihedral symmetry combined with unitarity im- 
ply a parameterization of scattering matrix in terms of 
three real variables. For two-dimensional square lattice, 
the most general form of S'-matrix is given by a 4x4 ma- 
trix: 




rA 

tAB 



tsA 


tcA 


toA \ 






tcB 


toB 




tsc 


rc 


tuc 




tBD 


tcD 


ro / 





(Al) 



where A, B, C, D are the coefficients of incoming plane 
waves. Primed values denote coefficients of outgoing 
waves. At this stage, one has 16 complex parameters 
for this S'-matrix. 

Unitarity condition {S^S = I) combined with time- 
inversion symmetry {S^^ = S*) gives — S (notice that 
5* is the transposed matrix, not the conjugate). It leaves 
10 complex parameters. Using four reflections of two 
types (1) A ^ B, and {2} A ^ C, B ^ D that generate 
the dihedral symmetry group D4 consequently reduces 
this number to three complex variables. We obtain the 
S'-matrix in the form (|37|l . Unitarity allows finally to 
express the scattering matrix with only 3 real parameters: 

p-H2|ti|2 + |t|||2^1 



rtr.+r*t\\ +2\t^\ 







Subtracting the third equation from the first one allows 
us to define a first real parameter ip: 



There remains two independent equations: 
I f2(|t|||2 + |t^|2) ^t||e-#+t*e'^ 




2\t» = 2Re[(i|| +t±)e- 



2|i|| -tj_\^ = 2Re[(i|| -ii)e-''^] 



{m,n+ll 

A. a 



{in,n+l} 



B 



{m,n+l ) 



A 



m+l,n) 



B, 



A 



+ 

B 



B 



m+l,n) 



tb': 



f 

{m,n-l} 



FIG. 7: In 2D case, each node is indexed by a pair of numbers 
{m,n}. Incoming and outgoing plane waves are connected by 
the 4x4 scattering matrix. 



Two more real parameters are needed to complete the 
parametrization: 

ft,! = cos0„e*(^"+'^) 
+t±= cos0„e*('^"+'^) 

Expressions of transmission and reflection coefficients as 
functions of these three real parameters are given in the 
main text, see Eq. (|38|l . We remark an interesting fact: in 
the case of perfect transmission {r — 0), only the separate 
thread solution (|t|| \ = l,t± = Q) is possible. 



APPENDIX B: BAND STRUCTURE FOR A 
SQUARE LATTICE OF WIRES 

In this appendix we derive the band-structure for a 
square lattice of wires of non-interacting electrons. As in 
the one dimensional case, the wave-function away from 
impurities (i.e. nodes here) could be written as combina- 
tion of plane waves: 



b: 



-ikx 



where i — {x, y} and the 

are defined 
By definition of the scattering matrix: 



(Bl) 

coefficients 
on Fig. 13 
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/im+l/ 



4 m,n+l 



= 5 



Dm.n+1 




Bloch periodicity condition for the wave-function implies: 



An.m 
An.rn 



^i(k'^—k)na^ik'^ma j^.Q 
^i{k' —k)m,a ik' na /lO,0 

e y e ^ Ay 

k'^na pO,0 

y 



ma^ik'yi 



Using the last two expressions, we obtain the secular equation: 



1 



1 



^^gi{k-k'^)a ^^^i{k'y-k'^+2k)a 



i(k — k'y)a 

r 



1 



„2ika 



^^i(k'^+k)a _ ^ 



= 



(B2) 



(B3) 



(B4) 



Replacing scattering matrix elements by their parametrization H38|l . we get the implicit band-structure equation given 
in the main text (|41(l . Here we propose two more different ways to write the same dispersion relation, where the 
combinations ka + if), k'^a and k'yU have been replaced respectively by the simpler notations fc, k'^ and fc^: 



cos(A; + 0„) cos(fe + (l)y) + cosfc^ cosk'y cos0„ cos^i, — - (cos /c^ + cos fc^)(cos(^„ cos(fc + (j)v) + cos</)„ cos(fc + 0„)) 



C0S(A: + (/)„)— cos 0„ cos fc^ COs(fc + — COS01, cosfc^ 

cos(fc + (j)u) — COS 4>u COS fc(^ cos(A; + 0„) — cos 0„ cos fc^ 

The first form is useful to identify symmetries of band-structure. The second form is useful to derive RGF equations 
directly without using the formalism developed in section ITlll We remark that in the ID case {t± = 0), and in the 
2D symmetric case {t± = <||) the band-structure equations are the same, namely: cos(fca + tp + </)) = cos (/> cos{k' a) 



APPENDIX C: ANY LATTICE 
GENERALIZATION 

In this part we will discuss particular points met in 
section Uni of this article. First of all we could obtain the 
dispersion relation for any network i.e. when the wires 
lengths are not necessarily equal. In that case Eq. H23|) 
is modified into: 

A,, = ^''^ e^'^^-^/^^li^e^'^^-Z^A,™ (CI) 

rn 

We choose the origin of coordinates needed to define the 
amplitudes Aij at the centers of each link. This for- 
mula means that the amplitude of the wave going from 
node j to node i is the sum of amplitudes coming from 
all neighbors m of node j, multiplied by phase factors 
exp{ikLjm/2) due to propagation from the middle of link 
(jm) to the node j, then scattered on node j with prob- 

ability amplitude 5"^;^ and finally reaching the middle of 
link (ij) with a new phase factor exp(ifcLy/2). We will 



now write the same equation in vector form. The ex- 
pression will be more transparent and this permits us 
to express the secular equation for energy eigenvalues 
/c in a compact form. If we fix the energy of the sys- 
tem then the stationary states are completely determined 
by 2Nl amplitudes, where is the number of links. 
The factor 2 arises since each wave can propagate in two 
opposite directions on each link. So the set of ampli- 
tudes {Ay} could be presented as a vector \ A) in a 2Nl- 
dimensional Hilbert space. We choose the orthonormal 
basis associated with network links {mn\ij) = SmiSnj- 
Each link is represented by two basis vector \ij) and 
this orientation difference should be taken into account 
in various summations over first neighbors. We define 
the vector \A) = ^ Aij\ij) and the length operator 

<ij> 

L = The vector form of Eq. IjCip reads: 

<ij> 



\A(k)) = e'''^/^fe'''^/^\A{k)) (C2) 
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The possible values of k are given by: 



det(e-*'^'^ - T) = 



(C3) 



One remarks that the periodicity of the spectrum 1)26(1 
is lost in the general case, unless there exists a Afc such 
that exp(iAfcL) = 1. Let us remind that this periodicity 
of the spectrum allowed us to evaluate the integral I{k) 
in Eq H29(l : the contribution of each band was the same 
and we replaced the sum over any filled band just by sum 
over all the eigen-vectors of operator T. 

The second point to be clarified is the spectral decom- 
position J2q 19)6"*'° (ql = T. It holds only if \q) vectors 



form an orthonormal basis. Orthogonality is clear as \q) 
is an eigen- vector of an unitary operator. 

(q\f\k) ^ e-*«"(q|fc) = e-''"'{q\k) (C4) 

Let us now evaluate the vector norm in the basis: 

{q\q)= E('^I*^')(*^19)= E (C5) 



<ij> 



<ij> 



But we know that the norm of wavefunction H22(l in the 
physical Hilbert space should be equal to unity. 



\ip{x)\'^dx = 1 = a 



network 



E 

<ij> 



\A,,{q)\'+A,,{q)A*M 



r 



sin(fca) 
ka 



(C6) 



So if we demand {q\q) — 1 we are doing an approximation 
neglecting the term proportional to sin(fca)/fca. This ap- 
proximation is legitimate in our case, as we consider sys- 
tems where the typical number ef electrons along each 
link between two nodes is large. Clearly, it will break 
down for links of the order of the Fermi wave-length. 
Supposing that {q\q) = 1 is equivalent to identify the 
norm in the (infinite dimensional) physical Hilbert space, 
with the norm associated to the orthonormal basis \ij) in 
the {2Nl dimensional) auxiliary Hilbert space. The fact 



that our equations are not depending explicitly on the 
network scale parameter a is closely related to this ap- 
proximation. So if one were to estimate finite size correc- 
tions to the RGF equation, one should take the physical 
normalization of the |g)-basis into account. Such cor- 
rections would likely produce RGF equations where the 
nodes on the lattice are no longer renormalized indepen- 
dently of each other, by contrast to what we obtained in 
section lim see Eq. 
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